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Q-[ Abstract 

It is proposed that all flavor mixing is caused by the mixing of the three quark and lcpton 
families with vectorlike fermions in 5 + 5 multiplcts of SU (5). This simple assumption implies that 
both Vckm and Umns are generated by a single matrix. The entire 3x3 complex mass matrix 
of the neutrinos M u is then found to have a simple expression in terms of two complex parameters 
and an overall scale. Thus, all the presently unknown neutrino parameters are predicted. The best 
fits are for 9 a t m ~ 40°. The leptonic Dirac CP phase is found to be somewhat greater than n. 

It is a striking fact that the leptonic mixing angles of the MNS matrix [1] are much larger than 
the corresponding quark mixing angles of the CKM matrix [2]. Grand unification suggests a simple 
explanation for this. In SU(5), a family of quarks and leptons is contained in the multiplets 10 + 5, 
with the left-handed leptons contained in the 5 and the left-handed quarks contained in the 10. Thus, 
if there is more mixing among the 5 multiplets of different families than among the 10 multiplets, the 
disparity between leptonic and quark mixing angles would be explained. This idea can be implemented 
in models based on any grand unified group, since all such groups contain SU (5) as a subgroup. Several 
ways of implementing this basic idea have been proposed in the literature [31 0]. 

Here we propose a model in which the three 10 + 5 families of fermions are supplemented by three 
5 + 5 pairs. (The possible existence of such additional "vectorlike" fermions has been much discussed 
in the literature in a variety of contexts [3 EJ El El IB]-) The central idea of the model proposed here is 
that all inter-family mixing is caused by the mixing between the 5 multiplets of the ordinary families 
and the 5 multiplets of the additional vectorlike pairs. As a consequence of having this common 
source, both quark mixing and lepton mixing are controlled in this model by a single matrix, which 
we call A. This matrix can be determined from the masses and mixing angles of the quarks alone, and 
this allows the entire 3x3 complex mass matrix M v of the known neutrinos (which contains 9 real 
physical observables) to be predicted in terms of just two complex parameters and an overall mass 
scale. The resulting formula turns out to be quite simple. In the "flavor basis" of the neutrinos, i.e. 
the basis (f e , u^, u T ), the neutrino mass matrix is given by 
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where pe ia , qe t/3 , and the overall scale \i v are free parameters of the model, S is the Kobayashi- 
Maskawa CP phase, and the Vy are the CKM matrix elements. The expression V us stands for 
s\n9 us cos6 us = V US \J\ — |V^ S | 2 . (Since the other quark mixing angles are small, their cosines can be 
set to one, and we can use simply V c b and V u b-) Note that the CKM elements Vij describing quark 
mixing enter this formula in such a way that if they vanish the neutrino mass matrix becomes diagonal 
and the neutrino mixing vanishes also. But while these CKM mixing parameters are small, they are 
here multiplied by large ratios of quark masses, explaining naturally why the neutrino mixing is of 
order one. (^-|F US | ~ 4, ^|V c6 | ~ 2, and '^\V ub \ ~ 3.) We shall see later how this arises. 

In this model, there are three families of fermions denoted by 10j + 5j (i = 1,2,3), and three 
vectorlike pairs of fermion multiplets denoted by 5^+5^ (A = 1, 2, 3). In the absence of the vectorlike 
pairs, the Yukawa couplings and mass matrices of the three families would be flavor diagonal, due to 
discrete symmetries that distinguish the three families from each other. All flavor mixing is indirectly 
caused by mass terms that mix the 5j with the 5 A . The model is defined by the following quark and 
lepton Yukawa terms: 

C Yu k = ri(10il0i)(5 H ) + y i (10 i 5 i )(5^) 

+ y i (10 i 10 i )(45 ff ) + y i (10 i 5 i )(45 t ff ) 



(2) 

(K.K.\/K„\/K„\ 

Mr' 



+ Y AB (5> A 5 B )(1 H ) + 24(5^X1^), 

where the subscript H denotes Higgs multiplets. The renormalizable Yukawa terms in Eq. (2) are the 
most general allowed by the symmetry K\ x K2 x K% x K', where (for a given i equal to 1,2, or 3) K{ is 
a Z2 symmetry under which 10«, 5j, and l' Hi are odd and all other fields even. K' is a Zn symmetry 
(N > 2) under which 5' A -> e 27ri / N 5' A , t A -> e 2ni / N t A , 1 H -> e - 47ri / N l H , and l' m -> e- 2ni / N l' Hi . 
(This is just one example. Many other simple discrete symmetries would give the form in Eq. (2).) 

The first four terms in Eq. (2) are the standard Yukawa terms of SU(5) grand unification, and are 
the minimal terms needed to give mass to the known quarks and leptons. (As already noted in the 
original Georgi-Glashow paper on SU(5) unification, the presence of a 45 of Higgs fields avoids the 
unrealistic relations between down quark and charged lepton masses that would arise if only a 5 of 
Higgs fields existed [9].) The fifth term is the standard dimension-5 Weinberg operator that gives the 
left-handed neutrinos Majorana masses |10| . (The symmetry K' prevents other dimension-5 operators 
that would give neutrino masses, such as 5 a 5 b 5h5h-) Note that all these standard terms are forced 
to be flavor diagonal by the Ki symmetries. 

The last two terms in Eq. (2) are the only ones peculiar to this model. The first of these simply 
gives masses to the vectorlike fermions, and the second gives masses that mix these vector like fermions 
with the three families. The Higgs fields in these two terms are gauge singlets, so that their VEVs 
would naturally be superlarge. All that matters for the purposes of this paper is that the masses 
coming from these two terms be roughly of the same scale, which we shall call A4, and that this scale 
be large compared to the masses of the down quarks and charged leptons. Note that the Yukawa 
matrices in these two terms are in general not diagonal. 

This is the model; all that remains is to extract its predictions. First, let us examine the mass 
matrix of the down quarks that emerges from Eq. (2). There are left-handed anti-down quarks in 
both 5j and 5 A , which will be denoted d\ and D A respectively. There are left-handed down quarks 
in both 10j and 5 A , which will be denoted di and D' A respectively. Altogether, then, there is a 6 x 6 
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mass matrix for the down quarks, given by 

<W> " « ( ( 7f m" b ) ( % ) ■ < 3 > 

where (mo); = 2/i(5jy) + ^(45^), Mas = Y^lff), and A^j = y' Aj (l' Hj ). Here and throughout the 
paper, Dirac mass matrices are multiplied from the left by the left-handed fermions and from the right 
by the right-handed fermions (or, equivalently, the left-handed anti-fermions) 

The 6x6 matrix in Eq. (3) can be block-diagonalized by multiplying it from the right by a unitary 
matrix U whose elements are of order one (since the elements of the matrices A and M are of the 
same order) and from the left by a unitary matrix whose angles are of order nip/M. <C 1 and which 
therefore can be neglected. Specifically, the unitary matrix U is such that 

m D \ f m D \ ( A B \ _ { m D A m D B \ 

A M ) U = { A M ) \ C D ) ~ \ M J' (4j 

where A = [I + A^M -1 tM" 1 A]- 1 /2. The off-diagonal block rri£)B in the last matrix in Eq. (4) can 
be removed by a rotation from the left that is of order nip/M, which is negligible, as already noted. 
Thus, after block-diagonalization, the upper-left 3x3 block that describes the masses of the three 
observed down quarks becomes simply 



M D = m D A. (5) 

In other words, the net effect of the mixing of the three families with the heavy vectorlike fermions 
is to multiply the diagonal mass matrix ran from the right by a non-diagonal matrix A. This can be 
understood diagramatically from Fig. 1. 
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Fig. 1 Diagrams showing how the mass terms (A) that mix the 5j with the 5^ lead to insertions of 
the matrix A on external 5 fermion lines. 

From these diagrams, it is easy to see that such factors of A accompany external fermion lines that 
are in the 5 representation of SU(5). Thus the mass matrices for the up quarks, down quarks, 
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charged leptons and neutrinos (which come, respectively, from (10 10)5//, (10 5)5#, (5 10)5-^, and 
(5 5)5#5.h- terms) have the form 



M v = mu, 



M D = m D A, 
M e = A T m e , 
M v = A T m v A 



(6) 



where the matrices mu, mp, me, and m„ are all diagonal. This can also easily be shown by block- 
diagonalizing the full mass matrices of the charged leptons and neutrinos in the same way that we did 
for the down quarks. 

It might seem that the matrix we have called A should be different for the different types of 
fermions due to renormalization effects. At the unification scale, the same 3x3 matrices A and M 
appear in the 6x6 mass matrices of the charged leptons and the down quarks. But, due to gluon 
loops, the A and M of the down quarks should run more strongly between the unification scale and 
the scale A4 than the corresponding matrices of the leptons. The crucial point, however, is that A 
depends on the ratio M _1 A; and since gauge boson loops cause A and M to run in the same way, 
these effects cancel out in A. Moreover, the renormalization effects due to Yukawa couplings (which 
are small for the 5 fermions) can be neglected. Thus, it really is the same matrix A that appears in 
Md, Me, and M v . Ultimately, this is due to SU(5) symmetry. It is this fact that makes this model 
so predictive. 

To extract the predictions of this model, let us consider first the quarks. The up quarks are already 
in the "mass basis" , since Mjj = mu is diagonal. The down quark mass matrix can be written 
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where fj, d is the 33 element of the diagonal matrix mo . As we have indicated in the second line of Eq. 
(7), the matrix A can be brought to triangular form by rotations of the right-handed quarks. Then, 
by rescaling the parameters [i d , e d , and 6 d , one can make the diagonal elements of A equal to 1. And 
finally, by rephasings of the left-handed and right-handed quarks, one can remove all phases except 
one, which we can place in the 13 element, as shown. 

Since the up quark mass matrix is already diagonal, the CKM mixing matrix of the quarks comes 
entirely from diagonalizing the mass matrix in Eq. (7). It is easy to show from that equation that 
\V cb \ = sm9 cb = tanflcfe = e d a, \V ub \ = sm9 ub = tan6> Mfc = 5 d c, \V US \^1 - \V US \ 2 = tan# us = 6 d b/e d , 
and 9 = 5, the CP phase of the quarks. (Notice that the basic structure predicted by the model, given 
in Eq.(6), explains why \V ub \ ~ |KsK*|, since \V US \ = 0(5 d /e d ), \V cb \ = O(e^), and \V ub \ = 0(5 d ).) 

It is also clear from Eq. (7) that e d = m s /m b and S d = m d /m b . Notice that this implies that S d <C 
< 1, so that the matrix mo is not only diagonal (because of family symmetry) but "hierarchical". 
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(Though the model as presented in Eq. (2) does not explain the mass hierarchy among the families, 
we shall see that a simple extension of the model can do this.) Combining the above relations gives 
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Turning to the leptons, one sees that the mass matrix of the charged leptons can be written 
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We have transformed A to have the same form as in the second line of Eq. (7), by doing the same 
combination of rotations to the left-handed leptons as we did to the right-handed down quarks, followed 
by analogous rescalings and rephasings. If we do the same rotations to the left-handed charged leptons 
and left-handed neutrinos, no MNS mixing is induced at this stage; but by doing so, it is clear that we 
also make the matrix A appearing in M v = A T m u A have the same form as in the second line of Eq. 
(7) and Eq. (9). What results is precisely the form shown in Eq. (1) for the neutrino mass matrix. 
Since we only have the freedom to rephase the left-handed neutrinos, there are three physical phases 
in Eq. (1), rather than one as in the other mass matrices. The extra two phases are the ones called 
a and (3 in Eq. (1). 

The hierarchy among the charged lepton masses tells us that 1, So, the diagonal 

matrix mi is hierarchical, just as mr> and m\j are. By comparing Eqs. (7) and (9), we see how this 
model explains the disparity between the neutrino mixing angles and quark mixing angles. Because 
Md = mo A , whereas Mi = A T mi, the mass matrix of the down quarks has a hierarchy among the 
rows, whereas the charged lepton mass matrix has a hierarchy among the columns. Since rotations 
of the left-handed fermions (which are the ones relevant to the CKM and MNS mixing angles) are 
rotations among the rows, we see that small quark mixing angles and large lepton mixing angles arise. 
(This is a realization of the basic idea of "lopsided" models OH].) 

The charged lepton mass matrix in the second line of Eq. (9) is not yet diagonal, but to a very 
good approximation it can be diagonalized by rotations done only to the right-handed charged leptons. 
Rotations of the left-handed charged leptons are also required, but they are rotations by angles that 
are proportional to e|, (5i/ei) 2 , and <5^. The only one of these that is numerically significant is a 



rotation in the \xlt~l plane by angle 6 



= ae 2 e ~ 



2(m At /m r ) 



0.4°. This contributes at the 1% level 



to the atmospheric neutrino mixing angle, and this has only a minor effect on the predictions of the 
model, as we shall see. 

Thus, in effect, the mass matrix in the second line of Eq. (9) is in the mass basis of the left-handed 
charged leptons. Consequently, the mass matrix M v shown in Eq. (1) contains all the information 
about the masses, mixings and CP-violating phases of the neutrinos. There are nine physical observ- 
ables involved: the three neutrino masses, the three MNS angles, the Dirac CP phase, and the two 
major ana CP phases of the neutrinos. These are all determined by the five model parameters in Eq. 
(1): /Ujy, pe ia , and qe l/3 . 



Since five neutrino observables have already been measured (9 so i, 6 a tm, #13, <5"?-23 an d Sw^), we 
can use them to determine the five model parameters, and then predict the four as-yet-unmeasured 
neutrino observables. Since the equations are non-linear (they involve trigonometric functions), there is 
no guarantee that they can fit the five measured neutrino properties with five adjustable parameters. 
(To put it another way, the adjustable parameters a and /3 are bound within the range [0, 
Nevertheless, good fits are obtained. This is only true, however, if some of the measured neutrino 
properties have values that lie within a smaller range than that presently allowed by experiment. 
For example, although the current experimental range of the atmospheric neutrino mixing angle is 
Gatm = 45 ± 6.5° [TT) . the model only obtains good fits for 6 atm ~ 43°, with values near 40° preferred, 
as we shall see. The fits also prefer a value of 9 so i greater than or equal to 34°, i.e. greater than 
the present experimental central value. The quark properties are also constrained: the best fits are 
obtained with m s /m,i ~ 20, and 5 greater than or equal to its present experimental central value. 
Thus, in addition to predicting the four as-yet-unmeasured neutrino observables, the model places 
non-trivial and testable constraints on the values of quantities that have been measured. 

In Table I, we show a representative fit in which all the input quark parameters and the neutrino 
observables obtained as output are in their experimentally allowed ranges (and in most cases at their 
central values). The experimental values are taken from the 2012 Review of Particle Properties |11| . 
except for Si ep (the neutrino Dirac CP phase) where we use the result of a recent global analysis of 
neutrino data [12]. For mb/m s we have used the renormalization group results of |13| to obtain m s (m b ) 
from m s (2GeV), which is given in |11| . 

Table I. A fit to the quark and neutrino data. ii u: pe ia , and qe % ^ are model parameters. 6i ep is the 
neutrino Dirac CP phase, and (M u ) ee the mass that comes into neutrinoless double beta decay. 
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2.32 x 10~ 3 eV 2 


2.32^ x 10~ 3 eV 2 




7.603 x 10~ 5 eV 2 


(7.5 ± 0.2) x 10~ b eV 2 


hep 


1.15-7T rad 


l.lnl^Z ^d 


{M v ) ee 


0.0020 eV 





Note that the model's prediction for 5[ ep is 1.157T radians, which accords remarkably well with the one- 
sigma range found in JT2J of 1.17t1q |^ rad. The value of (M u ) ee (to which the amplitude of neutrinoless 
double beta decay is proportional) is much smaller than the experimental limits, which tend to be in 
the range of a few tenths of an eV to several eVs for different experiments [11] . This prediction of the 
model is not very sensitive to variation of the model's input parameters. 

Figs. 2-4 show the degree of sensitivity of the 5i ep prediction to the values of 9 a tm, &soh an d S 



(the quark CP phase). In Fig. 2, we have fixed the values of all the quark mass ratios and CKM 
parameters, and of 9 so i and #13, but have allowed 9 a tm and the ratio Sm^/Sm^ (which we henceforth 
call r) to take different values. The curves are the relation of r to the predicted 5[ ep for different 
values of 9 a tm- The horizontal lines are the one-sigma limits for r. One sees that Si ep is predicted to 
be roughly 1.157T radians and that values of Q a tm ~ 41° are preferred. 

In Fig. 3, we have done a similar thing, but this time fixing 9 a tm to be 40° and allowing 6 so i and 
r to vary. One can see a preference for values of 6 so i equal or above the present experimental central 
value. In Fig. 4, we have allowed the quark CP phase 5 and r to vary. One sees that the best-fit value 
of 5[ ep is rather insensitive to the assumed values of the measured quark and neutrino properties, but 
the width of the range of 5i ep values that give good fits is quite sensitive. A more detailed analysis of 
the predictions of the model will be given in another paper. 
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Fig. 2 The result of fits with the values of quark parameters given in Table I, 6 so i = 34.2°, and 
#13 = 9.12°. The curves are the relation of r (= Sm^/Sm^) to the predicted 5i ep for different values 
of 

9atm- The horizontal lines are the one-sigma limits for r. 
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Fig. 3 The result of fits with the values of quark parameters given in Table I, 9 a tm = 40°, and 
#13 = 9.12°. The curves are the relation of r (= 5m^ 2 / 'Sm^) to the predicted 8\ ep for different values 
of 9 soi . 
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Fig. 4 The result of fits with the values of quark parameters other than 5 given in Table I, 9 so i = 34.2°, 
and 6»i3 = 9.12°, and 

@atm — 40°. The curves are the relation of r (= bm^l '^^23) ^° the predicted 
5i ep for different values of the quark phase 5. 

Besides its great simplicity, one feature of the model proposed here that increases its plausibility 
is that it allows the simple solution to the Strong CP Problem proposed in [6]. All that is required is 
that CP be assumed to be a symmetry of the Lagrangian that is spontaneously broken by the VEVS 
of the singlet scalars l' Hi that produce the off-diagonal mass matrix A in Eq. (3). In fact, the model 
proposed here has the same structure as the model originally proposed by A.E. Nelson in [6], except 
that here the 3x3 mass matrices of the ordinary three families are required to be flavor diagonal by 
family symmetries. 

The model proposed here gives an account of how the CKM and MNS flavor mixings arise, but 
does not explain the mass hierarchy among the families, since the hierarchies in the diagonal matrices 
mcr, uid, rri£, and m v are simply assumed. There are, however, several simple ways in which the 
present model could be extended to give an explanation of the mass hierarchy. One way is to combine 
the structure in this model with the structure assumed in [3]. In that paper, the mass hierarchies 
were explained by the three ordinary families mixing with vectorlike 10 + 10 fermion pairs in a way 
analogous to the mixing with 5 + 5 assumed here. Combining the structures of the two models would 
be appealing since it would mean that the vectorlike fermions would comprise entire family-antifamily 
pairs. (It has been pointed out that this can lead in a simple way to unification of gauge couplings in 
non-SUSY models [8].) 

Another possibility would be a Froggatt-Nielsen scheme [14 1. For example, instead of the family 
symmetry assumed above, suppose that the i*Q (i = 1 or 2) were Z4 symmetries, under which 5j — > — 5j, 
l' Hi — > — l' Hi , 10i — > ilOj, and S% — > iSi, where the Si are Froggatt-Nielsen fields that are SU(5) 
singlets. Effective Yukawa terms containing factors of 10j would then have to contain equal numbers 
of factors of {Si) /Mp = e^, where Mp is some flavor-physics scale. If ei <C €2 <C 1, a mass hierarchy 
among families would result. Moreover, the hierarchy would be strongest for the up quarks (for which 
it is quadratic in the e^), intermediate for the down quarks and charged leptons (for which it is linear 
in the e^), and weakest for the neutrinos (which involve no factors of e,). This is just the pattern that 
is observed. One should note that the same relationship among the mass hierarchies is obtained in 
the approach of [3] . 

Finally it should be noted that the present model could be embedded in many grand unified 
schemes. For example, in an 5O(10) model, the ordinary families could be in three 16 multiplets, 
while the vectorlike fermions could be in three 10 multiplets. In E$, one gets the extra vectorlike 
fermions "for free", since the 27 contains 16 + 10 + 1 of 5*0(10). Different patterns of breaking of 
the grand unified group could be assumed without affecting the predictions for fermion masses and 




mixings. For example, in many unified models, an adjoint Higgs field does some of the breaking of 
the unification group. If that adjoint Higgs multiplet does not transform under the K' symmetry 
mentioned after Eq. (2), it would not couple renormalizably to {h' A h' B ) or (5^5j) and hence not 
contribute to the matrices A and M in Eq. (4) and the matrix A. Consequently, except for negligible 
higher-order corrections, the matrix A would not "know" that the unification group is broken, and 
the same A would appear in both the quark and lepton sectors, as is necessary for the model to be 
predictive. 

In conclusion, if all flavor changing in both the quark and lepton sectors arises as a consequence 
of the mixing of the ordinary families with vectorlike fermions that are in 5 + 5 of SU(5), a testable 
relationship arises between the quark and lepton mixing. This relationship allows the prediction of 
the four as-yet-unmeasured neutrino observables as well as testable constraints on several quantities 
that have been measured. Measurement of the Dirac CP phase of the neutrinos 5i ep , as well as more 
precise determinations of such quantities as 6 atm , so ii \V u b\, m s/ m d, and 6 (the quark CP phase) 
would provide stringent tests of the model. 
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